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Abstract 



We consider the quantum complexity of estimating matrix elements of unitary irreducible represen- 
tations of groups. For several finite groups including the symmetric group, quantum Fourier transforms 
yield efficient solutions to this problem. Furthermore, quantum Schur transforms yield efficient solu- 
' tions for certain irreducible representations of the unitary group. Beyond this, we obtain poly(n)-time 

quantum algorithms for approximating matrix elements from all the irreducible representations of the al- 
ternating group A„, and all the irreducible representations of polynomial highest weight of U{n), SU{n), 
and SO{n). These quantum algorithms offer exponential speedup in worst case complexity over the 
fastest known classical algorithms. On the other hand, we show that average case instances are classi- 
cally easy, and that the techniques analyzed here do not offer a speedup over classical computation for 
the estimation of group characters. 



1 Introduction 



I Explicit representations of groups have many uses in physics, chemistry, and mathematics. All representa- 

. tions of finite groups and compact Unear groups can be expressed as unitary matrices given an appropriate 

\^ I choice of basis^. This makes them natural candidates for implementation using quantum circuits. Here we 

. show that polynomial size quantum circuits can implement: 

o 

• The irreducible representations of any finite group which has an efficient quantum Fourier transform. 
This includes the symmetric group Sn- 

• The irreducible representations of the alternating group An- 

• The irreducible representations of polynomial highest weight of the unitary U{n), special unitary 
■ SU{n), and special orthogonal SO{n) groups. 

^ . Using these quantum circuits one can find a polynomially precise additive approximation to any matrix 

element of these representations by repeating a simple measurement called the Hadamard test, as described 
in section [2l 

More precisely, for the finite groups Sn and An we obtain any matrix element of any irreducible rep- 
resentation to within ±e in time that scales polynomially in 1/e and n. For the Lie groups U{n), SU{n), 
and SO{n) we obtain any matrix element of any irreducible representation of polynomial highest weight 
to within ±e in time that scales polynomially in 1/e and n. Because the representations considered are of 
exponentially large dimension, one cannot efficiently find these matrix elements by classically multiplying 
the matrices representing a set of generators. Note that, many computer science applications use multiplica- 
tive approximations. In this case, one computes an estimate i of a quantity x with the requirement that 
(1 — e)x < X < (1 — e)^;. The approximations obtained in this paper are all additive rather than multiplicative. 



'Parts of this work were completed at MIT's Center for Theoretical Physics and RIKEN's Digital Materials Laboratory. 
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Figure 1: The complexity results on the symmetric group refer arbitrary irreducible representations in Young's 
orthogonal form. The results on the braid group refer to the Jones- Wenzl representations, which give rise to Jones 
and HOMFLY polynomials. The complexity class DQCl is the set of problems solvable in polynomial time on a one 
clean qubit computer. It is generally believed that one clean qubit computers are weaker than standard quantum 
computers but still capable of solving some problems outside of BPP. 



For some problems, the computational complexity of additive approximations can differ greatly from that of 
mulitplicative approximations [TT], [2] . 

For exponentially large unitary matrices, the typical matrix element is exponentially small. Thus for 
average instances, a polynomially precise additive approximation provides almost no information. However, 
it is common that the worst case instances of a problem are hard whereas the average case instances are 
trivial. In section [5] I narrow down a class of potentially hard instances for the problem of additively 
approximating the matrix elements of the irreducible representations of the symmetric group to polynomial 
precision. I also present a classical randomized algorithm to estimate normalized characters of the symmetric 
group Sn to within ±e in poly(n, 1/e) time. (The character is normalized by dividing by the dimension of 
the representation, so that the character of the identity element of the group is 1.) Thus, the techniques 
described here for evaluating matrix elements of irreducible representations of groups on quantum computers 
do not provide an obvious quantum speedup for the evaluation of the characters of 5„. 

Our results on the symmetric group relate closely to the quantum complexity of evaluating Jones polyno- 
mials and other topological invariants. Certain problems of approximating Jones and HOMFLY polynomials 
can be reduced to the approximation of matrix elements or characters of the Jones- Wenzl representation 
of the braid group, which is a g-deformation of certain irreducible representations of the symmetric group 
[3l [37l [33l [25] . Figure [T] compares the complexity of estimating matrix elements and characters of the Jones- 
Wenzl representation of the braid group to the complexity of the corresponding problems for the symmetric 
group. Exact complexity characterizations {i.e. completeness results) are not known for all of these problems, 
and the exact relationships between the complexity classes referenced in figure [T] are not rigorously known. 
Nevertheless, the results seem to suggest that in general the matrix elements are harder to approximate 
than the normalized characters, and that the Jones- Wenzl representation of braid group is computationally 
harder than the corresponding irreducible representations of the symmetric group. 



2 Hadamard Test 

The Hadamard test is a standard technique in quantum computation for approximating matrix elements of 
unitary transformations. Suppose we have an efficient quantum circuit implementing a unitary transforma- 
tion U, and an efficient procedure for preparing the state | ■(/')• We can then approximate the real part of 
{tP\U lip) using the following quantum circuit. 



H — 



The probability of measuring |0) is 



Po = 



1+Remum 



Thus, one can obtain the real part of (V'l U jV') to precision e by making 0(l/e^) measurements and counting 
what fraction of the measurement outcomes are |0). Similarly, if the control bit is instead initialized to 
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-^(|0) — i |1)), one can estimate the imaginary part of (-01 U Thus the problem of estimating matrix 
elements of unitary representations of groups reduces to the problem of implementing these representations 
with efficient quantum circuits. 

3 Fourier Transforms 

Let G be a finite group and let G be the set of all irreducible representations of G. We choose a basis for 
the representations such that for any p £ G and g G G, 5 is represented by a dp x dp unitary matrix with 
entries pi_j{g). The quantum Fourier transform over G is the following unitary operator [29j 

dp rj~ 

Here \g) is a computational basis state (bitstring) indexing the element g of G. Similarly, \p,i,j) is three 
bitstrings, one indexing the element p € G, and two writing out the numbers i and j in binary. The standard 
discrete Fourier transform is the special case where G is a cyclic group. 
The regular representation of any g € G is 

U, = J2 \9h) {h\ . 

heG 

A short calculation shows 

dp dp 
peG »J = 1 i' J' = l 

In other words, by conjugating the regular representation of g with the quantum Fourier transform, one 
recovers the direct sum of all irreducible representations of g~^. 

Given an efficient quantum circuit implementing Uft one can thus efficiently estimate any matrix element 
of any irreducible representation of G using the Hadamard test. Quantum circuits implementing the Fourier 
transform in polylog(|G|) time are known for the symmetric group 8 and several other groups [28\ The 
matrix elements of the representations depend on a choice of basis. The bases used in quantum Fourier 
transforms are subgroup adapted (see [28] ). In particular, the symmetric group Fourier transform described 
in [H] uses the Young- Yamanouchi basis, also known as Young's orthogonal form. 

In section [S] we describe a more direct quantum circuit implementation of the irreducible representations 
of the symmetric group, which generalizes to yield efficient implementations for the alternating group. 

4 Schur Transform 

Let H be the Hilbert space of n d-dimensional qudits. 

H = (C*)^". 

We can act on this Hilbert space by choosing an element u G U{d) and applying it to each qudit. 

We can also act on this Hilbert space by choosing an element tt G S'„ and correspondingly permuting the n 
qudits. 
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M®" and Af^ are reducible unitary nd-dimcnsional representations of U{d) and Sn, respectively. These two 
actions on Ti. commute. 

The irreducible representations of Sn are in bijective correspondence with the partitions of n. Any 
partition of n into d parts indexes a unique irreducible representation of U{d). U{d) has infinitely many 
irreducible representations, so these partitions only index a special subset of them. As discussed in fB^, there 
exists a unitary change of basis C/schm- such that 

A 

where A ranges over all partitions of n into d parts. 

As shown in [B], J7schur can be implemented by a poly(n, d) size quantum circuit. Thus, using the 
Hadamard test, one can efficiently obtain matrix elements of these representations of the symmetric and 
unitary groups. 

5 Complexity of Symmetric Group Representations 

As described in section [31 quantum computers can solve the following problem with probability 1 — (5 in 
poly(n, 1/e, log(l/i5)) time. Note that standard Young tableaux index the Young- Yamanouchi basis vectors, 
as discussed in section ISTTI 

Problem 1: Approximate a matrix element in the Young- Yamanouchi basis of an irreducible representation 
for the symmetric group Sn- 

Input: A Young diagram specifying the irreducible representation, a permutation from S'„, a pair of 
standard Young tableaux indicating the desired matrix element, and a polynomially small parameter e. 
Output: The specified matrix element to within ±e. 

It appears that no polynomial time classical algorithm for this problem is known. Due mainly to applica- 
tions in quantum chemistry, many exponential time classical algorithms for the exact computation of entire 
matrices from representations of the symmetric group have been developed [551 UHl [3H1 [3S1 [SI [HI [311 150] . 
There appears to be no literature on the computation or approximation of individual matrix elements of 
representations of Sn- 

On the other hand, the precision of approximation achieved by the quantum algorithm is trivial for average 
instances. We can see this as follows. Let A be a Young diagram of n boxes, let px be the corresponding 
irreducible representation of S'„, and let d\ be the dimension of p\. For any n d Sn, the root mean square 
of the matrix elements of px{n) is 

y ^ a.beB 

where B is any complete orthonormal basis for the vector space on which p\ acts. We see that 

Y,\{a\pxin)\b)\'^l 

aeB 

since, by the unitarity of px{7r), this is just the norm of \b). Thus, 



RMS (pa W) = 




The interesting instances of problem 1 are those in which d\ is exponentially large. In these instances, 
the typical matrix element is exponentially small, by equation [1] Running the quantum algorithm yields 
polynomial precision, thus one could instead simply guess zero every time, with similar results. 
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Figure 2: Here is a sequence of Young diagrams, such that as the number of boxes increases, the Young diagram 
converges asymptotically to some fixed shape, in this case a triangle. 



That the average case instances are trivial does not mean that the algorithm is trivial. Hard problems 
that are trivial on average are a common occurrence. The most relevant example of this is the problem 
of estimating a knot invariant called the Jones polynomial. A certain problem of estimating the Jones 
polynomial of knots is BQP-complete [TBI [S] [1]. The Jones polynomial algorithm is based on estimating 
matrix elements of certain representations of the braid group to polynomial precision. On average these 
matrix elements are exponentially small. Nevertheless, the BQP-hardness of the Jones polynomial problem 
shows that the worst-case instances are as hard as any problem in BQP. 

By analogy to the results on Jones polynomials, one might ask ask whether problem 1 is BQP-hard. 
The existing proofs of BQP-hardness of Jones polynomial estimation rely on the fact that the relevant 
representations of the braid group are dense in the corresponding unitary group. Thus, one can construct 
a braid whose representation implements approximately the same unitary as any given quantum circuit. 
Furthermore, it turns out that the number of crossings needed to achieve a good approximation scales only 
polynomially with the number of quantum gates in the circuit. Unlike the braid group, the symmetric 
group is finite. Thus, no representation of it can be dense in a continuous group. Hence, if the problem of 
estimating matrix elements of the symmetric group is BQP-hard, the proof will have to proceed along very 
different lines than the BQP-hardness proof for Jones polynomials. 

Lacking a hardness proof, the next best thing is to identify a class of instances in which the matrix 
elements are large enough to make the approximation nontrivial. As shown below, we can do this using 
the asymptotic character theory of the symmetric group. Note that we need not worry about the matrix 
elements being too large, because even if we know a priori that a given matrix element has magnitude 1, it 
could still be nontrivial to compute its sign. 

Let TT be a permutation in Sn, and let A be a Young diagram of n boxes. The character 

XA(7r) = Tr(pA(7r)) 

is clearly independent of the basis in which px is expressed. Furthermore, the character of a group element 
depends only on the conjugacy class of the group element, because for any representation p, 

T^{p{hgh-^)) = T^{p{h)p{g)p{hr^) = Tr(p(g)). 

To understand the behavior of the characters of Sn as n becomes large, consider a sequence of Young 
diagrams Ai, A2, A3, . . ., where A„ has n boxes. Suppose that the diagram A„, when scaled down by a factor 
of l/y/n, converges to a fixed shape cu in the limit of large n, as illustrated in figure [2l Let c?a„ be the 
dimension of the irreducible representation corresponding to Young diagram A„. Let tt be a permutation in 
Sk- We can also consider tt to be an element of Sn for any n > k which leaves the remaining n ~ k objects 
fixed. As shown by Biane^SJ, 

^^f^ = aHn-l-l/2-|-0(n-l^l/2-i). (2) 

Here |7r| denotes the minimum number of transpositions needed to obtain tt. Note that these are general 
transpositions, not transpositions of neighbors. C^(w) is a constant that only depends on tt € Sk and the 
shape uj. A precise definition of what it means for the sequence to converge to a fixed shape is given in [9], 
but for present purposes, the intuitive picture of figure [5] should be sufficient. 
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Xa„ (tt) / dx^ is the average of the matrix elements on the diagonal of (tt) . In the present setting, where 
TT is fixed, Xa„ (tt) /(iA„ shrinks only polynomially with n. Thus polynomial precision is sufficient to provide 
nontrivial estimates of these matrix elements. Nevertheless, finding diagonal matrix elements of (tt) for 
fixed TT and large n is not computationally hard. This is because, as discussed in section [51 the Young- 
Yamanouchi basis is subgroup adapted. Thus, for any tt which leaves all bit the first k objects fixed, px„ (tt) 
is a direct sum of irreducible representations of tt in 5*^. Because k is fixed, any irreducible representations of 
Sk has dimension 0(1) and can therefore be computed in 0(1) time by multiplying the matrices representing 
transpositions. 

To produce a candidate class of hard instances of problem 1, we recall that the character xx^ (ti") depends 
only on the conjugacy class of vr. Thus, we consider tt' conjugate to tt. Like tt S Sn, tt' G 5„ leaves at least 
n — k objects fixed, and the representations XA„(7r') have diagonal matrix elements with polynomially small 
average value. However, the objects left fixed by tt' need not be k + 1, k + 2, . . . ,n. Indeed, tt' can be chosen 
so that the object n is not left fixed, in which case px^iji') cannot be written as the direct sum of irreducible 
representations of Sm for any m < n. 

There is an additional simple way in which an instance of problem 1 can fail to be hard. Let r(7r) be the 
minimal number of transpositions of neighbors needed to construct the permutation tt. If r(7r) is constant 
or logarithmic, then the matrix elements of the irreducible representations of tt can be computed classically 
in polynomial time by direct recursive application of equation 1131 For a class of hard instances of problem 
1 I propose the following. 

Hypothesis 1 Let t: be a permutation in Sn- We consider it to permute a series of objects numbered 
1, 2, 3, . . . , n. Let s(7r) be the number of objects that ir does not leave fixed. Let /(tt) be the largest numbered 
object that tt does not leave fixed. Let r(7r) be the minimum number of transpositions of neighbors needed to 
construct tt. Let X be a Young diagram of n boxes, and let px be the corresponding dx- dimensional irreducible 
representation of Sn- I propose the problems of estimating the diagonal matrix elements of px{tt) such that 
s(7r) = 0(1), '(tt) — i^{n), and r(7r) = fl{n) as a possible class of instances of problem 1 not solvable 
classically in polynomial time. 

Although this hypothesis contains many restrictions on vr, it is clear that permutations satisfying all of these 
conditions exist. One simple example is the permutation that transposes 1 with n. 

6 Characters of the Symmetric Group 

Because characters do not depend on a choice of basis, the computational complexity of estimating characters 
is especially interesting. Hepler 24] showed that computing the characters of the symmetric group exactly is 
#P-hard. It is clear that an algorithm for efficiently approximating matrix elements of a representation can 
aid in approximating the corresponding character. Specifically, the quantum algorithm for problem 1 yields 
an efficient solution for the following problem. 

Problem 2: Approximate a character for the symmetric group 5*7^. 

Input: A Young diagram A specifying the irreducible representation, a permutation vr from Sm and a 
polynomially small parameter e. 

Output: Let x^{t^) be the character, and let dx be the dimension of the irreducible representation. The 
output Xout must satisfy |xout — X^i"^) I dx\ < e with high probability. 

However, as we show in this section, problem 2 is efficiently solvable using only classical randomized 
computation. Thus the techniques used for problem 1 do not offer immediate benefit for problem 2. Although 
this is in some sense a negative result, it provides an interesting illustration of the difference in complexity 
between estimating individual matrix elements of representations and estimating the characters. 

We can reduce problem 2 to problem 1 by sampling uniformly at random from the standard Young 
tableaux compatible with Young diagram A. For each Young tableau sampled we estimate the corresponding 
diagonal matrix element of pxi'^)^ as described in problem 1. By averaging the diagonal matrix elements 



6 



for polynomially many samples, we obtain the normalized character to polynomial precision. The problem 
of sampling uniformly at random from the standard Young tableaux of a given shape is nontrivial but it 
has been solved. Greene, Nijenhuis, and Wilf proved in 1979 that their "hook-walk" algorithm produces the 
standard Young tableaux of any given shape with uniform probability [H]. Examination of [H] shows that 
the time needed by the hook-walk algorithm to produce a random standard Young tableaux compatible with 
a Young diagram of n boxes is upper bounded by O(n^). 

By averaging over diagonal matrix elements we lose some information contained in the individual matrix 
elements. This observation gives the intuition that it should often be harder to estimate individual matrix 
elements of a representation than to estimate its trace. Jones polynomials provide an example in which this 
intuition is confirmed. As discussed in [33], computing the Jones polynomial of the trace closure of a braid 
reduces to computing the normalized character of a certain representation of the braid group. The problem 
of additively approximating this normalized character is only DQCl-complete. In contrast, the individual 
matrix elements of this representation yield the Jones polynomial of the plat closure of the braid and are 
BQP-complete to approximate. We see a very similar phenomenon in the symmetric group; problem 2 is is 
solvable by a randomized polynomial-time classical algorithm, whereas problem 1 is not, as far as we know. 

To construct a classical algorithm for problem 2, first recall that the character of a given group element 
depends only on the element's conjugacy class. We can think of any tt e 5'„ as acting on the set {1, 2, . . . , n}. 
The sizes of the orbits of the elements of {1, 2, . . . , n} under repeated application of tt form a partition of 
the integer n. For example, consider the permutation tt G ^5 defined by 

7r(l) = 2 7r(2) = 3 7r(3) = 1 7r(4) = 5 7r(5) = 4. 

This divides the set {1,2,3,4,5} into the orbits {1,2,3} and {4,5}. Thus it corresponds to the partition 
(3, 2) of the integer 5. Two permutations in Sn are conjugate if and only if they correspond to the same 
partition. Thus, we can introduce the following notation. For any two partitions and A of n define to 
be the irreducible character of 5„ corresponding to the Young diagram of A evaluated at the conjugacy class 
corresponding to fi. 

To obtain an efficient classical solution to problem 2 we use the following theorem due to Roichman[32j. 

Theorem 1 (From |32] ) For any partitions fi ~ (/ii, . . . , /i;) and A = (Ai, . . . , Afc) of n, the corresponding 
irreducible character of Sn is given by 

A 

where the sum is over all standard Young tableaux A of shape A and 

w^iA)^ n 

l<i<k 

where B{fj,) = {/ii + . ■ ■ + Hr\^ < r < 1} and 

{—1 box i + 1 of A is in the southwest of box i 
i + 1 is in the northeast of i, i + 2 is in the southwest of i + 1, and i + 1 ^ B^fjb) 
1 otherwise 

By using the hook walk algorithm we can sample uniformly at random from the standard Young tableaux 
A of shape A. By inspection of theorem [T] we see that for each A sampled we can compute ty^(A) classically 
in poly(n) time. By averaging the values of Wp(A) obtained during the course of the sampling we can 
thus obtain a polynomially accurate additive approximation the the normalized character, thereby solving 
problem 2. 

Some readers may notice that theorem[T]is similar in form to the much older and better-known Murnaghan- 
Nakayama rule. However, the Murnaghan-Nakayama rule is based on a sum over all "rim-hook tableaux" 
of shape A (see ,32:). It is not obvious how to sample uniformly at random from the rim-hook tableaux of 
a given shape. Thus, it is not obvious how to use the Murnaghan-Nakayame rule to obtain a probabilistic 
classical algorithm for problem 2. 
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7 Lie Groups 



7.1 Introduction 

Because U{n), SU{n) and SO{n) are compact linear groups, all of their representations are unitary given the 
right choice of basis pj. In section [4] we described how to efficiently approximate the matrix elements from 
certain unitary irreducible representation of U(n). Here we present a more direct approach to this problem, 
which can handle a larger set of representations of U(n) and also extends to some other compact Lie groups: 
SU{n) and SO{n). 

U{n), SU{n), and SO{n) are subgroups of GL{n), the group of all invertible n x n matrices. All of the 
irreducible representations of U{n) and SU (n) can be obtained by restricting the irreducible representations 
of GL{n) to these subgroups. The best classical algorithms for computing irreducible representations of 
GL{n) and U{n) appear to be those of [H] and [50] . These classical algorithms work by manipulating matrices 
whose dimension equals the dimension of the representation. Thus, they do not provide a polynomial time 
algorithm for computing matrix elements from representations whose dimension is exponentially large. The 
implementation of irreducible representations of 5*0(3) and SU{2) by quantum circuits has been studied 
previously by Zalka 40 . 

7.2 Gel'fand-Tsetlin representation of U{n) 

The irreducible representations of the Lie group U{n) are most easily described in terms of the corresponding 
Lie algebra u{n). It is not necessary here delve into the theory of Lie groups and Lie algebras, but those who 
are interested can see [TH]. For now it suffices to say that u{n) is the set of all antihermitian nx n matrices, 
and for any u G U{n) there exists h S u{n) such that u = e''. Given any representation a : uin) u{m) 
one can construct a representation A : Uin) U{m) as follows. For any u € U{rL) find a corresponding 
h{u) £ u{n) such that e'^ — u, and set A{u) — e''^'''^"^^. If a is an antihermitian representation of u{n) 
then A is a unitary representation of U{n). Furthermore, it is clear that A is irreducible if and only if a is 
irreducible. 

It turns out that the irreducible representations of the algebra glin) of all n x n complex matrices remain 
irreducible when restricted to the subalgebra u{n). Furthermore, all of the irreducible representations of 
uin) are obtained this way. Let Eij be the n x n matrix with all matrix elements equal to zero except for 
the matrix element in row i, column j, which is equal to one. The set of all n^ such matrices forms a basis 
over C for gl{n). Thus to describe a representation of gl{n) it suffices to describe its action on each of the 
Eij matrices. 

As described in chapter 18, volume 3 of [35] . explicit matrix representations of gl{n) were constructed by 
Gel'fand and Tsetlin. (See also [H].) In their construction, one thinks of the representation as acting on the 
formal span of a set of combinatorial objects called Gel'fand patterns. The Gel'fand-Tsetlin representations 
of Ep^p-i and i?p-i,p are sparse and simple to compute for all p e {2,3,..., n}. This property makes the 
Gel'fand-Tsetlin representations particularly useful for quantum computation. 

A Gel'fand pattern of width n consists of n rows of integer^ The j'^ row (from bottom) has j entries 
"^ij j"^2j , ■ ■ • I'Ti'j.j- (Note that, in contrast to matrix elements, the subscripts on the entries of Gel'fand 
patterns conventionally indicate column first, then row.) These entries must satisfy 

mj^n+l > "^j\n ^ 'rrij+i.n+i- 

Gel'fand patterns are often written out diagrammatically. For example the Gel'fand pattern of width 3 with 
rows 

mi^3 = 4 TO2,3 = 1 m3,3 = 

TOi^2 = 3 TO2,2 = 
?71l,l = 2 

1 Some sources omit the top row, as it is left unchanged by the action of the representation. 
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is represented by the diagram 



This notation has the advantage that the entries that appear directly to the upper left and upper right of a 
given entry form the upper and lower bounds on the values that entry is allowed to take. 

We call the top row of a Gel'fand pattern its weighi^. To each weight of width n corresponds one ir- 
reducible representation of gl{n). This irreducible representation acts on the formal span of all Gel'fand 
patterns with that weight (of which there are always finitely many) . To describe the action of the represen- 
tation of gl{n) on these patterns let 



•'P,q 



bi 





nr=i(^j,p 




-i) nr=i ('i,p-2 - ij,p-i - 1) 


1 ^ 


I\i^j{h,p- 


-1 — 


^j;P~l){h,p~l ~ ^j,p-l ~ 1) 






^3:P- 


-1 + 1) I\i=li^i,p-2 - lj,p-l) 


1 = 


Ui^j{h,p- 


-1 — 


1 ) 1 1 ~^ 



1/2 



1/2 



(3) 
(4) 

(5) 



Let M be a Gel'fand pattern and let M^^ be the Gel'fand pattern obtained from M by replacing nij^p with 
TT^j,p "I" 1. Similarly, let Mp be the Gel'fand pattern in which Trij p has been replaced with irij p — 1. The 
representation a,fi of gl(n) corresponding to weight m G Z" is defined by the following rule^, known as the 
Gel'fand-Tsetlin formulas. 



(6) 



am{Ep-i,p)M 


j=i 


afi,{Ep,p^i)M 


= E^P-i^P-i 
j=i 


am{Ep^p)M 


/ p p~i 

= E"^-p~E 

j=i 



M 



(7) 



(8) 



These formulas give implicitly a representation for all of gl{n), because any Eij can be obtained from 
operators of the form Ep^i p and by using the commutation relation [Eik,Eki] = Eu. By restricting 

the representation a,fi to antihermitian subalgebra of gl{n) and taking the exponential, one obtains an 
irreducible group representation : U(n) — > U{d^), where is the number of Gel'fand patterns with 
weight rh. 

It should be noted that some references claim that the set of allowed weights for representations of 
GL(n) is N", whereas others identify, as we do, Z" as the allowed set of weights. The reason for this 
is that irreducible representations of GL{n) in which the entries r7i„^i,TO„^2j . ■ . ,n^n,n of the weight are all 
nonnegative are polynomial invariants [26]. That is, for any g G GL{n) and any m € N", each matrix element 
of the representation Pffi{u) is a polynomial function of the "n? matrix elements of u. The representations 
involving negative weights are called holomorphic representations, and many sources choose to neglect them. 
In the case that m S N", the Gel'fand diagrams of width n bijectively correspond to the semistandard Young 
tableaux of n rows (c/. [15], pg. 517). 

■^It is actually the highest weight of the representationfSS]. but for brevity I just call it the weight throughout this paper. 
^Warning: |35l contains a misprint, in which the sums in equations [6l and Ffl are taken up to j = p instead of jr = p — 1. 
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7.3 Quantum Algorithm for U(n) 

In this section we obtain an efficient quantum circuit implementation of any irreducible representation of 
U{n) in which the entries mi^„, . . . ,mn^n of the highest weight are all at most polynomially large. The 
dimension of such representations can grow exponentially with n. Unlike the Schur transform, the method 
here does not require TOi,n, ■ • ■ , to be nonnegative. We start by finding a quantum circuit implementing 
the Gel'fand-Tsetlin representation of an n x n unitary matrix of the form 



uo = 



where all off-diagonal matrix elements not shown are zero. After that we describe how to extend the 
construction to arbitrary n x n unitaries. 

For a given weight m € Z" we wish to implement the corresponding representation Ajj^{uq) with a 
quantum circuit. To do this, wc first find an n x n Hcrmitian matrix Hq such that e'^° = uq. It is not hard 
to see that Hq can be computed in polynomial time and takes the form 

/ill hi2 

^22 



_ 

Thus, 

Hq = hiiEii + /112-E12 + ^12^21 + /122-E22- (9) 

Hence, 

am{Ho) = hiia^-,{Eii) + /ll2am(£^12) + ^t2am(£^2l) + h22a^{E22)- (10) 

To implement 74,j;(uo) with a quantum circuit, wc think of a.ffi{HQ) as a Hamiltonian and simulate the 
corresponding unitary time evolution e"""^*^^"^* for t = —1. The Hamiltonian a^{Hfj) has exponentially 
large dimension in the cases of computational interest. However, examination of equation [9] shows that Hq 
is a linear combination of operators of the form Ep^p^i and i?p_i_p. Thus, by the Gel'fand-Tsctlin rules 
of section [7.21 affi{HQ) is sparse and that its individual matrix elements are easy to compute. Under this 
circumstance, one can use the general method for simulating sparse Hamiltonians proposed in [4 . 

Define row-sparse Hamiltonians to be those in which each row has at most polynomially many nonzero 
entries. Further, define row-computable Hamiltonians to be those such that there exists a polynomial time 
algorithm which, given an index i, outputs a list of the nonzero matrix elements in row i and their locations. 
Clearly, all row computable Hamiltonians are row-sparse. As shown in the unitary e"'^* induced by any 
row-computable Hamiltonian can be simulated in polynomial time provided that the spectral norm \\H\\ and 
the time / are at most polynomially large. We have already noted that a^{Ho) is row-computable. a^{Ho) 
is row sparse, and because we are considering only polynomial highest weight, the entries of the Gcl'fand 
patterns, and hence the matrix elements of g^^Hq) are only polynomially large. Thus, by Gershgorin's circle 
theorem ||ara(i/o)|| is at most poly(n). 

Having shown that a quantum circuit of poly(n) gates can implement the Gel'fand-Tsetlin representation 
of an n X n unitary of the form uq, the remaining task is to extend this to arbitrary n x n unitaries. 
Examination of the preceding construction shows that it works just the same for any unitary of the form 

Up — Ip © li © llj^_p_2, 



Ull U12 
"21 "22 

1 

1 
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where Ip denotes the p x p identity matrix and it is a 2 x 2 unitary. Corresponding to Up is again an 
antihermitian matrix of the form 



Hp ^ Op® hi 



where Op is the pxp matrix of all zeros and ft, is a 2 x 2 antihermitian matrix such that — u. The only issue 
to worry about is whether |ja,7i(i?p)|| is at most poly(n). By symmetry, one expects that ||a,7i(7fp)|| should 
be independent of p. However, this is not obvious from examination of equations [3] through [8] Nevertheless, 
it is true, as shown in appendix [Xj Thus, the norm is no different than in the p — Q case, i.e. Hq. 

By concatenating the quantum circuits implementing Aj^(ui), A^{u2), ■ ■ • one can implement 

A^{uiU2 . . . Ml). We next show that any n x n unitary can be obtained as a product of poly(n) matrices, 
each of the form Up, thus showing that the quantum algorithm is completely general and always runs in 
polynomial time. 

For any 2x2 matrix M, let £{M,i,j) be the n x n matrix in which M acts on the i*"^ and j*"^ basis 
vectors. In other words, the k, I matrix element of £{M, i,j) is 



Mil 

Afi2 

M22 
5kt 



if 
if 

if A: 
if A: 



i and I — i 
i and / — j 
j and / = i 
j and I = j 



otherwise 



Thus 



Next note that. 



= £ 



( 


Ull 


"12 




"21 


"22 



, m + 1 , m + 2 



"11 "12 
"21 "22 

1 

1 



m + 1, m + 3 I = 
, TO + 2, m + 3 ) 5 



Thus the matrix 



"11 "12 

"21 "22 

"11 "12 

"21 "22 



is obtained as a product of three matrices of the form Up. By repeating this conjugation process, one can 
obtain 

"11 "12 
"21 "22 



£ 



,TO+ 1, 171 + 2 £ 



, TO + 1, TO + 3 



1 

1 



, m + 2, TO + 3 



(11) 



for arbitrary i,j as a product of one matrix of the form 



( 


"11 


"12 




"21 


"22 



for some p and at most 0{n) matrices of the form 



£ 





■ 


1 " 


( 


1 






,P+l,P + 2 



,9+1,9 + 2 



with various q. A matrix of the form shown in equation [TT] is called a two-level unitary. As shown in section 
4.5.1 of [29 , any n x n unitary is obtainable as a product of poly(n) two-level unitaries. Thus we obtain 
A^(U) for any n x n unitary U using poly(n) quantum gates. One can then obtain any matrix element of 
AmiU) to precision ±e by repeating the Hadamard test 0(l/e^) times. 
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7.4 Special Orthogonal Group 

The special orthogonal group SO{n) consists of all n x n real orthogonal matrices with determinant equal 
to one. The irreducible representations of SO{n) are closely related to those of U{n) and can also be 
expressed unitarily using a Gel'fand-Tsetlin basis. As discussed in chapter 18, volume 3 of [3S], the nature 
of the representations of SO{n) depends on whether n is even or odd. Following 35] and |T5], we therefore 
introduce an integer k and consider SO{2k + 1) and SO{2k) separately. 

The irreducible representations of SO{2k + 1) are in bijective correspondence with the set of allowed 
weight vectors m consisting of k entries, each of which is an integer or half-integer. Furthermore, the entries 
must satisfy 

mi^n > m-2,n > ■ • • > mk^n > 0. 

The irreducible representations of SO{2k) correspond to the weight vectors m with fc — 1 entries, each of 
which must be an integer or half integer, and which must satisfy 

As in the case of U{n), the set of allowed Gel'fand patterns is determined by rules for how a row can 
compare to the one above it. For SO{n) these rules are slightly more complicated, and the rule for the j"^ 
row depends on whether j is odd or even. Specifically the even rule for j = 2k is 



™l,2fe+l > ™l,2fc > ™2,2fc+l > "i2,3fc > ■ ■ ■ > Wfe,2fc+1 > "ifc^a/c > -™fc,2fc-l, 

and the odd rule for j — 2k ~ 1 is 

m-l,2fc > "T,i_2fc-1 > m2,2fc > '7l2,2fc-l > • • • > ™fc-l,2fc > mk-l,2k~l > |™fc,2fc|- 

The Lie algebra so{n) corresponding to the Lie group SO{n) is the algebra of all antisymmetric n x n 
matrices. For any G € SO(n) there exists a 5 e so(n) such that = G. The Lie algebra so{n) is the space 
of all n X n real traceless antisymmetric matrices. Thus it is spanned by operators of the form 



IkA — E; 



'i,k 



k,i 



Ki <k < 



We can fully specify a representation of so{n) by specifying the representations of the operators of the form 
Iq+i,q because these generate so{n). That is, any element of so{n) can be obtained as a linear combination 
of commutators of such operators. The Gel'fand-Tsetlin representation 6^ of these operators depends on 
whether q is even or odd, and is given by the following formulas. 



G2p{M) = 
bfh{l2p+i,2p)M = 



rir^l [(^r,2p-l - \Y ~ ik 



lj.2p 



r,2p+l 



1)2 
2' 



J.2p 



r#jAS-,2p 
r,2p 



2p 



i,2p+ 



-P MP 

1) nr=l(^r,2p+2 ^ ^|,2p+l) 



1/2 



P 



l(4'j,2p+l 



i)n. 



p..^,){p 



r^j\'-r,2p+l ''j,2p+lJ\''j,2p+l 



2p+l 



1/2 



n 



p+i 7 

1 tr,2p-|-2 



nr=l ^r,2p-|-l(^r,2p+l ^ 1) 



p 



-nM^p' 



2p+2,2p+l 



]M 



P 



2p+ 



i(M)M+; 



E^2p+l(^^2p^l)^^2p^l 



iC2p{M)M 



By applying these rules to the set of allowed Gel'fand patterns described above one obtains the irreducible 
representations of the algebra so{n). By exponentiating these, one then obtains the irreducible represen- 
tations of the group SO{n). Thus the quantum algorithm for approximating the matrix elements of the 
irreducible representations of SO{n) is analogous to that for U{n). 
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7.5 Special Unitary Group 



The irreducible representations of SU{n) can be easily constructed from the irreducible representations of 
U{n), using the following facts taken from chapter 10 of [7]. The representations of U{n) can be partitioned 
into a set of equivalence classes of projcctively equivalent representations. Two representations of U{n) 
with weights I — (/i, ^2, ■ ■ • , In) E^nd m — {mi, m2, . . . , are projcctively equivalent if and only if there 
exists some integer s such that rrii = h + s for all 1 < i < n. Any irreducible representation of U{n) 
remains irreducible when restricted to SU{n). Furthermore, by choosing one representative from each class 
of projcctively equivalent representations of U{n) and restricting to SU{n) one obtains a complete set of 
inequivalent irreducible representations of SU{n). The Lie algebra su{n) corresponding to the Lie group 
SU{n) is easily characterized; it is the space of all traceless n x n antihermitian matrices. Thus the matrix 
elements of the irreducible representations of SU{n) are obtained by essentially the same quantum algorithm 
given for U{n) in section [7751 



7.6 Characters of Lie Groups 

As always, an algorithm for approximating matrix elements immediately gives us an algorithm for approx- 
imating the normalized characters. However, the characters of U{n), SU{n), and SO(n) are classically 
computable in poly(n) time. As discussed in |17j . the characters of any compact Lie group are given by the 
Weyl character formula. In general this formula may involve sums of exponentially many terms. However, 
in the special cases of U{n), SU{n), and SO{n) the formula reduces to simpler forms [T7]. given below. 

Because characters depend only on conjugacy class, the character Xmiu) depends only on the eigenvalues 
of u. For u £ Uin) let Ai, . . . , A„ denote the eigenvalues. Let rfi — (mi,TO2, . . . ,?7i„) € Z" be the weight of 
a representation of U (n) . Let 

li = rrii + n — i (12) 
for each i G {1,2, . . . ,n}. The character of the representation of weight m is 



detA 
detB 



where A and B arc the following n x n matrices 



B,, = K-'- 

This formula breaks down if u has a degenerate spectrum. However, the value of the character for 
degenerate u can be obtained by taking the limit as some eigenvalues converge to the same value. As 
shown in [36] , one can obtain the dimension of the representation corresponding to a given weight m by 
calculating lim„^ix^(u). Specifically, by choosing Xj — e*-'' for each 1 < j < n and taking the limit as 
e — > one obtains 



where U is as defined in equation 1121 

As discussed in section ITTSl the irreducible representations of SU{n) are restrictions of irreducible rep- 
resentations of U{n), therefore the characters of SU{n) are given by the same formula as the characters of 
U{n). 

SO{n) consists of real matrices. The characteristic polynomials of these matrices have real coefficients, 
and thus their roots come in complex conjugate pairs. Thus, the eigenvalues of an element g € SO{2k -\- 1) 
take the form 

Ai, A2, . . . , Afc, 1, A^^, A2, . . . , A^, 
and for g e SO{2k), the eigenvalues take the form 

Ai, A2, . . . , Afc, Aj, A2, . . . , AJl. 
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As discussed in [T7] , the characters of the special orthogonal group are given by 

so(2fe+i), . ^ detC 
det D 

and 

so(2k), . ^ det g + det 
~ det G 

where C and D are the following fc x fc matrices 

^ -.m^+n— z+l/2 \ —(mj +n— 2+1/2) 

^ yi-i+l/2 _ ^-(ri-,:+l/2) 

and -E, F, G are the following (fc — 1) x (fc — 1) matrices 

where is as defined in equation [T^ 

As with U(n), the character of any element with a degenerate spectrum can be obtained by taking an 
appropriate limit. 



7.7 Open Problems Regarding Lie groups 

The quantum circuits presented in the preceeding sections efficiently implement the irreducible representa- 
tions of U{n), SU{n), and SO{n) that have polynomial highest weight and polynomial n. It is an interesting 
open problem to implement irreducible representations with quantum circuits that scale polynomially in the 
number of digits used to specify the highest weight. Alternatively, one could try to implement an Schur 
transform to handle exponential highest weight, which is also an open problem. It is even concievable that 
Schur-like transforms could be efficiently implemented for exponential n. That is, there could exist a quan- 
tum circuit of polylog(n) gates implementing a unitary transform V such that for any U G U(ri), VUV~^ is 
a direct sum of irreducible representations of U . Of course, if n is exponentially large, than we cannot have 
an explicit description of U , rather the group element U could itself be defined by a quantum circuit. 

A completely different open problem is presented by the symplectic group. Having constructed quantum 
circuits for SO{n) and SU{n), the symplectic group is the only "classical" Lie group remaining to be 
analyzed. Thus it is natural to ask whether its irreducible representations can be efficiently implemented by 
quantum circuits. Two different groups can go by the name symplectic group depending on the reference. 
Connected non-compact simple Lie groups have no nontrivial finite-dimensional unitary representations (see 
[7], theorem 8.1.2). This applies to one of the groups that goes by the name of symplectic. On the other 
hand, the irreducible representations of the compact symplectic group seem promising for implementation 
by quantum circuits. The main task seems to be finding a basis for these representations that is subgroup 
adapted and makes the representations unitary. A non- unitary subgroup-adapted basis is given in [27j . 



8 Alternating Group 

In section [31 we described a method to approximate matrix elements of the irreducible representations of 
the symmetric group using the symmetric group quantum Fourier transform. Here we take a more direct 
approach to this problem, which extends to the alternating group. To do this we must first explicitly describe 
the Young- Yamanouchi representation of the symmetric group. 
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Figure 3: The above matrices are irreducible representations in the Young- Yamanouchi basis with Young diagram 
EF°. Here Gi is the permutation in 54 that swaps i with i + 



8.1 Young- Yamanouchi Representation 

For a given Young diagram A, let Va be the vector space formally spanned by all standard Young tableaux 
compatible with A. For example, if 

A= 

then Va is the 3-dimensional space consisting of all formal linear combinations of 
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For any given Young diagram A, the corresponding irreducible representation in the Young- Yamanouchi 
basis is a homomorphism p\ from Sn to the group of orthogonal linear transformations on Va- It is not 
easy to directly compute p\{tt) for an arbitrary permutation tt. However, it is much easier to compute the 
representation of a transposition of neighbors. That is, we imagine the elements of Sn as permuting a set 
of objects 1, 2, . . . , n, arranged on a line. A neighbor transposition ai swaps objects i and j -I- 1. It is well 
known that the set {cti, (J2, . . . , Cn-i} generates Sn- 

The matrix elements for the Young- Yamanouchi representation of transpositions of neighbors can be 
obtained using a single simple rule: Let A be any standard Young tableau compatible with Young diagram 
A then 

Pa('^0A=^A+^1-^A', (13) 

where A' is the Young tableau obtained from A by swapping boxes i and and r/^ is the axial distance 

from box i -|- 1 to box i. That is, we are allowed to hop vertically or horizontally to nearest neighbors, and t 
is the number of hops needed to get from box i -I- 1 to box i, where going down or left counts as -f 1 hop and 
going up or right counts as —1 hop. To illustrate the use of equation 1131 some examples are given in figure 

El 

In certain cases, starting with a standard Young tableau and swapping boxes i and i + I does not yield 
a standard Young tableau, as illustrated below. 

standard nonstandard standard nonstandard 



i 


i+1 




i+1 


i 



i+1 i 
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Figure 4: To obtain the conjugate A of Young diagram A, reflect A about its diagonal. In other words the number of 
boxes in the i"^ column of A is equal to the number of boxes in the i**^ row of A. 



Some thought shows that all such cases are of one of the two types shown above. In both of these types, the 
axial distance is ±1. By equation [T51 the coefficient on the invalid Young tableau is y^T^^^^^^j^ — 0. Thus 
the representation lies strictly within the space of standard Young tableaux. 

8.2 Direct Quantum Algorithm for Sn 

We can directly implement the irreducible representations of Sn by first decomposing the given permutation 
into a product of transposition of neighbors. The classical bubblesort algorithm achieves this efficiently. For 
any permutation in 5„, it yields a decomposition consisting of at most 0{n^) transpositions. As seen in the 
previous section, the Young- Yamanouchi representation of any transposition is a direct sum of 2 x 2 and 1x1 
blocks, and the matrix elements of these blocks are easy to compute. As shown in [4J, any unitary with these 
properties may be implemented by a quantum circuit with polynomially many gates. By concatenating at 
most 0{n^) such quantum circuits we obtain the representation of any permutation in Sn- The Hadamard 
test allows a measurement to polynomial precision of the matrix elements of this representation. 

8.3 Algorithm for Alternating Group 

Any permutation tt corresponds to a permutation matrix with matrix element i,j given by The 
determinant of any permutation matrix is ±1, and is known as the sign of the permutation. The permutations 
of sign +1 are called even, and the permutations of sign —1 are called odd. This is because a transposition 
has determinant —1, and therefore any product of an odd number of transpositions is odd and any product 
of an even number of transpositions is even. 

The even permutations in Sn form a subgroup called the alternating group An, which has size n\/2. An 
is a simple group [i.e. it contains no normal subgroup) and it is the only normal subgroup of Sn other 
than {1} and Sn- As one might guess, the irreducible representations of the alternating group are closely 
related to the irreducible representations of the symmetric group. Consequently, as shown in this section, the 
quantum algorithm of section [521 can be easily adapted to approximate any matrix element of any irreducible 
representation of An to within ±e in poly(n, 1/e) time. 

Explicit orthogonal matrix representations of the alternating group are worked out in :34j and recounted 
nicely in [23j . Any representation p of Sn is automatically also a representation of A„. However an irreducible 
representation p of 5„ may no longer be irreducible when restricted to A„. Each irreducible representation 
of Sn either remains irreducible when restricted to An or decomposes into a direct sum of two irreducible 
representations of An. All of the irreducible representations of An are obtained in this way. 

The conjugate of Young diagram A is obtained by reflecting A about the main diagonal, as shown in figure 
|4l If A is not self-conjugate then the representation p\ of Sn remains irreducible when restricted to A„. In 
this case we can simply use the algorithm of section 18.21 If A is self-conjugate then the representation p\ 
of Sn becomes reducible when restricted to An. It is a direct sum of two irreducible representations of An, 
called PA+ and px-. The two corresponding invariant subspaces of the reducible representation are the -1-1 
and —1 eigenspaces, respectively, of the "associator" operator S defined as follows. 

Let A be a self-conjugate Young diagram of n boxes. Let Aq be the "typewriter-order" Young tableau 
obtained by numbering the boxes from left to right across the first row, then left to right across the second 
row, and so on, as illustrated in figure [S] For any standard Young tableau A of shape A, let S 5'„ be the 
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Figure 5: For a given Young diagram, there is a unique Young tableau in "typewriter" order, in which the boxes are 
numbered from left to right across the top row then from left to right across the next row, and so on, as illustrated 
in the example above. 



permutation that brings the boxes into typewriter order. That is, w\A = Aq. Let A be the conjugate of A, 
obtained by reflecting A about the main diagonal. If A is standard then so is A. Let d{\) be the length of 
the main diagonal of A. S is the linear operator on Va defined by 

An orthonormal basis for each of the eigenspaces of S can be easily constructed from the Young- 
Yamanouchi basis. When (n — c?(A))/2 is odd, every standard Young tableau A of shape A has the property 
sign(u;A) = — sign(u'^), and S" is a direct sum of 2 x 2 blocks of the form 

"0 -i ' 
i 

interchanging A and A. In this case, the linear combinations -^{A + iA) for each conjugate pair of standard 
Young tableaux form an orthonormal basis for the +1 eigenspace of S, and the linear combinations -^{A—iA) 
form an orthonormal basis for the —1 eigenspace of S. Similarly, when (n — d{\))/2 is even, sign(z«A) = 
sign(?«^) for all standard Young tableaux A of shape A. Thus S* is a direct sum of 2 x 2 blocks of the form 

"0 -1 " 
-1 

interchanging A and A. In this case the linear combinations -^(A — A) form an orthonormal basis for the 
+1 eigenspace of S and the linear combinations -^(A + A) form an orthonormal basis for the —1 eigenspace 
of 5. 

Suppose A is self-conjugate and {n — d{\)) /2 is even. Any matrix element of the irreducible representation 
p\+ of An is given by 

i(A + A)p,(7r)(r + f), 

where A, T is some pair of standard Young tableaux and tt is some element of A„. This is a linear combination 
of only four Young- Yamanouchi matrix elements of p\ (tt) . One can use the algorithm of section 18.21 to 
calculate each of these and then simply add them up with the appropriate coefflcients. The cases where 
in — d{X))/2 is odd and/or we want a matrix element of px- are analogous. 
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A ||a^(iL,)|| is independent of p 



As shown in section 17.31 the irreducible representation of an arbitrary u ^ U (n) with weight m can be 
computed by simulating the time evolution according to a series of Hamiltonians of the form A^{Hp), where 
is the Gel'fand-Tsetlin representation of the Lie algebra su{n) and 

Hp = Op® /i©0„_p_2, 

where /i is a 2 x 2 antihermitian matrix. The quantum algorithm for simulating these Hamiltonians require 
that ||A^(i?p)|j be at most poly(n). In section [731 we showed this to be the case for p — 0. Here we prove 
it for all p by showing: 

Proposition 1 Let h be a fixed 2x2 antihermitian matrix and let Hp — Op (B h (B On-p-2- Let a^i he the 
Gel'fand-Tsetlin representation of su(n) with weight rn. Then ||a/?i(ifp)|| is independent of p. 



Proof: 



Let = e^"^. Then 



e^^ ® l„-p-2- 



Thus for any < q < n, there exists V € U{n) such that 



(15) 



Specifically, V is just a permutation matrix. Let Affi be the Gel'fand-Tsetlin representation of SU{n). That 



Thus 



_d_ 

dfc 



= \MHp)e^-^^"-^\\ 
= \\affi{Hp)\\. 

Here we have used the fact that A^^ is a unitary representation. Similarly, 



(16) 



MH,)\\ 



d 



Using equation [T5l this is equal to 



^A^r^iVU^^V-') 



dfc 



Because A^-^ is a group homomorphism and V is independent of fc this is equal to 

ArniV) (^Aa^([/^^)^ A„-,{V)-' 

Because A^-j^ is a unitary representation this is equal to 

d 



dfc 



By equation [TBI this is equal to ||a^j(ffp)||. □ 
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